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ABSTRACT 

We report results of the (J = 3 —2) and HCO+ (J = 4 —3) observations 
of the W40 Hll region with the ASTE 10 m telescope (HPBW~ 22") to search 
for molecular outflows and dense clumps. We found that the velocity held in the 
region is highly complex, consisting of at least four distinct velocity components 
at WsR — 3, 5, 7, and 10 km s“^. The 7 km s ^ component represents the 
systemic velocity of cold gas surrounding the entire region, and causes heavy 
absorption in the spectra over the velocity range 6 < Wsr ^ 9 km s“^. 

The ~5 and ~10 km s“^ components exhibit high temperature (> 40 K) 
and are found mostly around the H ll region, suggesting that these components 
are likely to be tracing dense gas interacting with the expanding shell around the 
H II region. Based on the data, we identihed 13 regions of high velocity gas 
which we interpret as candidate outhow lobes. Using the HCO’*' data, we also 
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identified six clumps and estimated their physical parameters. On the basis of 
the ASTE data and near-infrared images from 2MASS, we present an updated 
three-dimensional model of this region. In order to investigate molecular outflows 
in W40, the SiO (J = 1 — 0, n = 0) emission line and some other emission lines 
at 40 GHz were also observed with the 45 m telescope at the Nobeyama Radio 
Observatory, but they were not detected at the present sensitivity. 

Subject headings: ISM: molecules-ISM:clouds-ISM: individual (W40) -stars: for¬ 
mation 


1. INTRODUCTION 


The W40 complex is a site of active star formation in the Aquila Rift containing molecu¬ 
lar clouds, star cluster, massive stars, and an H ll region. Several mid and far-infrared obser 
vations have revealed that the W 40 complex consists of a number of hlamentary structures 
to fo r m an hourglass shape fe.g.. iRodnev fc ReipurthI l2008l: iKuhn et all I2OIOI: I Andre et ah 


1^^ I— 

iects fYSOsl 

associated with this complex (e.g.. 

Rodriguez et al. 

2010: 

Kuhn et al. 

2010; 

Maurv et al. 

2011; 

Mallick et al. 

2013 

). The cluster and its associated Hll region are lo- 


cated to the northwest of the constriction of the hour glass. The cluster i s found to have a 
considerably large stellar surface density of 650 pc ^ flMallick et al.ll2013l) . Several ionizing 


sources fo r the H ll region are also found in this region. IRS lA South is an 09.5 spectral 


type star (IShuping et ahl 120121) and is the most luminous source. Another high mass star. 


IRS 5, is located in the northwest of IR S lA South. This is classihed as a BIV type star 


fjShuping et al.ll2012l: iMallick et al.ll201311 . Radi o emission from a rou ghly spherical compact 


Hll region has been detected around IRS 5 bv IMallick et ahl (120131) . They suggested that 
high-mass star formation is going on therein. 

Previous molecular line and radio co ntinuum observations revealed the presence of the 
dense gas around the W4 0 H ll region (e.g ., Zeilik fc Ladalll978 : Zhu et aPboOG : Pirogov et al. 
2013 : Maury et al. 2011 ). Dobashi et al. f 20051) and Dobashi ( 2011 ) named the dense region 
as TGU 279-P7 and No. 1195, based on the Digitized Sky Survey (DSS) and the Two Micron 
All Sky Survey (2MASS) extinction data, r espectively. The d e nse region has a ring -like 
shape along the boundary of the Hll region (IMaurv et al.ll201lh . iPirogov et ahl (l2013f) car¬ 
ried out multiple-molecular line observations such as HGO^, GS, and NH3, and found the 
large variations of the fractional abundances of these molecules over the dense region. All of 
these studies indicate that the H ll region is strongly interacting with the surrounding dense 
gas. The W40 complex is therefore an excellent region to study the interaction of massive 



































































































3 


stars with their parental molecular cloud. 


Several studies have discussed the geometry o f the W40 c omplex (e.g., Crutcherl 119771: 
Vallee 1987 : Kuhu et ah 2010 : Pirogov et ah 2013 ). Especially, Vallee ( 1987 ) suggested the 
blister model for W40 which is illustrated in their Figure 2 and is now regarded as a standard 
picture of the region. However, the complex velocity structure has prevented us from com¬ 
plete unde rstanding of the geometry and kine matics of the molecular cloud associated with 
W40 fe.g.. lCrutcherlll977l:IPirogov et al.ll2013[) . More detailed studies to uncover the internal 
density distribution and kinematics of the molecular cloud is needed to better understand 
the structure of this region. 


Low-mass members of this region were also revealed bv iMaurv et al.l (120111) and iMallick et al 


(120131) . They identihed a number of Class O/I sources and starl ess cores a r ound the Hll re¬ 
gion. However, outflow activity in this region remains uncertain. IZhu et al.l (120061) discovered 
a weak molecular outflow through their (J = 3 — 2) observations with the KOSMA 3 

m telescope whose half power beam width (HPBW) is ~ 80". The identihed outhow lobes 
appear to have complicated structure. In addition, the surveyed area is limited to a small 
region just around the H ll region. It is important to carry out a comprehensive survey for 
molecular outhows covering much wider area to reveal the star formation activity in this 
region. 


In this paper, we present results of the (J = 3 — 2) and HCO'*’ (J = 4 — 3) 

observations made toward a large area (~ 8' x 8') covering the whole W40 H ll region and its 
interacting molecular cloud, using Atacama Submillimeter Telescope Experiment (ASTE) 10 
m telescope. For the purpose of an outhow survey, we carried out (J = 3 — 2) mapping 
observations to detect 13 outhow candidates. In order to explore the velocity and spatial 
distributions of the dense gas in the W40 region, we also carried out HCO^ (J = 4 — 3) 
observations to identify six clumps, and further attempted to model their locations relative 
to the Hll region along the line-of-sight mainly through a comparison with the 2MASS 
near-infrared images. 


Hereafter we will refer to the observed region simply as W40. The distance estimation 
to W 40 remains uncertain, and varies from 340 to 600 pc depending on the methods used 


ng., Rodney fc Reipurth 2008 : Kuhn et al. 201ol: 


latest estimate of 500 pc derived by IShuping et al. 


Shuping et al.ll2012l) . Here, we adopt the 
()2012 ) throughout this paper. 
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2. OBSERVATIONS 


2.1. ^^CO and HCO'*^ Observations with the ASTE 10 m Telescope 


The mapping observations of W40 were performed with the (J = 3 — 2) and 

HCO'*' (J = 4 — 3) emission lin es at 345.795991 GHz and 356.73424 2 GHz, respectively, 
using the ASTE 10 m telescope flEzawa et ahll2004l: iKohno et ahll2004l) at Atacama, Chile. 
Data of these molecular lines were obtained simultaneously during the period between 2013 
September 26th and October 1st. The area observed is denoted in Figure [H which covers the 
central part of the W40 H ll region. The angular resolution of the telescope is HPBW~ 22" at 
the observed frequencies. A side-band separating (2SB) mixer receiver CATS345 was used for 
the observations. The spectrometer was an XF-type digital spectro-correlator, called MAC, 
with 1024 channels covering a bandwidth of 128 MHz. The frequency resolution was 125 
kHz, corresponding to a velocity resolution of 0.1 km s“^ at the observed frequencies. The 
system noise temperature Tgyg was 300 — 400 K during the observations. The observations 
were made using the on-the-fly (OTF) mapping tech nique fISawada etal.l 120081) . Data were 
calibrated with the standard chopper-wheel method f Kutner fc Ulichl 1981 ). The telescope 
pointing was regularly checked every 1.5 hour using ^^CO (J = 3 — 2) line from W-Aql, and 
the pointing error was smaller than 2" during the observations. In addition, we checked the 
intensity reprod ucibility by observ ing M17SW [R.A. (2000) = 18^20”^23^.08, Dec. (2000) 
=—16°11'43.3"] ( Wang et al. 1994) and Serp ens South [R.A. (2000) = 18^30^03®.8, Dec. 
(2000) =—02°03'03.9"] ( Nakamura et al. 2011) every day to calibrate the spectral data, which 
fluctuates within ~ 10%. 


Data reduction and map generation were performed using the software package NOS¬ 
TAR developed by the Nobeyama Radio Observatory (NRO). The data were convolved with 
a spheroidal function 0 and resampled onto a uniform grid of 10", and then corrected for 
the main-beam efficiency of the telescope (r^mb = 0.6). We finally obtained the spectral data 
in the brightness temperature scale Tmb in Kelvin with an effective angular resolution of 
approximately 31" and a velocity resolution of 0.1 km s“^ for each emission line. Typical 
noise level of the spectra was ATmb = 0.2 K for this velocity resolution. 


Schwab (1984) and Sawada et al. ( 20081) described the details of the spheroidal function. We applied the 


parameters, m = 6 and a = 1.0, which define the shape of the function. 
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2.2. SiO, CCS, HC 3 N, and HC 5 N Observations with the NRO 45 m Telescope 


We also carried out mapping observations toward W40 in the SiO (J = 1 — 0, u = 0) 
emission line at 43.42376 GHz using the NRO 45 m telescope on 2014 May 27th and 28th . 
The SiO molecule is a well-known probe of molecular outflows (e.g., iHirano et al.N2006l) . 
The molecule is also one of the shock tracers in molecular outflows which was hrst evidenced 


in LI 157 bv iMikami et al.l (119921). and has been wide l y used as the shock tr acer for other 
star forming regions (e.g., IShimaiiri et al.N2008l. 120091: iNakamura et al.l 120141) . In order to 
search for the possible SiO emission associated with the outflows in W40, we performed SiO 
mapping observations over the same regions observed with the ASTE telescope. 


We used a dual-polarization receiver Z45 flTokuda et al.ll2013l) and the digital spectrom¬ 
eters SAM45 with 4096 channels covering a bandwidth of 125 MHz. Combination of the 
receiver and spectrometers made it possible to observe the SiO, CCS (JAr =43 — 32 , 45.379033 
GHz), HC 3 N (J = 5-4, 45.490316 GHz), and HC 5 N (J = 17-16, 45.264721 GHz) emission 
lines simultaneously. The channel width of the spectrometers was 3.81 kHz, corresponding 
to a velocity resolution of ~ 0.025 km s“^ at 45 GHz. The observations were performed using 
the OTF mapping technique and the data were calibrated with the standard chopper-wheel 
method. Tgys was typically 180 K, and the pointing accuracy was better than 10" during 
the observations, as was checked by observing the SiO maser IRC-l-00363 at 43 GHz every 
1.5 hours, and HPBW were 0.7 and 37", respectively. Data reduction and analysis were 
done in the same manner as for the ASTE observations. We hnally obtained the spectral 
data with an effective angular resolution of approximately 57" and a velocity resolution of 
0.05 km s“^. The resulting 1 a noise level for this velocity resolution was AT^b ~ 0.2 K, 
but no signihcant SiO and other molecular emission was detected at this sensitivity. 


3. RESULTS AND DISCUSSION 
3.1. Distributions of Molecular Gas 


Figure [T] (a) shows the ^^CO (J = 3 — 2 ) intensity distribution integrated over the 
velocity range —3.0 A Ulsr = 15.0 km s“^ superposed on the DSS2-Red image. The intensity 
distribution in the south-west part of the map shows an arc- shaped structure which can be 


ident ihed also in other molecular maps in the literature (e.g., IZhu et al.ll2006l: iPirogov et al. 
2 OI 3 I) . Our map shows that the ^^GO emission is extending to the north-east direction. 


which has not been revealed by molecular observations with a sub-arcmin angular resolution 
to date. In the hgure, we also show the positions of known YSOs. i .e., C lass O/I sources (plus 
signs or circles) and starless cores (squares) found bv iMaurv et ahl (1201 if) through the 1.2 mm 







































6 


dust c ontinuum observations as well as Class I sources (circles) identified by iMallick et al. 
(120131) from the archival Spitzer observations in conjunction with near-infrared data from 
UKIRT. The high mass sources IRS lA South and IRS 5 are also shown with white stars in 
Figured] (a). It is clearly seen that the YSOs are distributed along the arc-shaped structure 
where the emission is enhanced up to Tmb(^^CO) = 40 — 50 K, suggesting that the 
arc-shaped structure is the dense gas heated and/or compressed by the ionization front of 
the H II region. 

Figure dl (b) shows the HCO’*' (J = 4 — 3) intensity distribution integrated over the 
velocity range 3.0 A Vlsr = 11-4 km s“^ superposed on the Herschel SPIRE 250 pm image 
B (colored in gray). The HCO'*’distribution is well correlated with t he dust distribution 


traced by the Herschel data and 1.2 mm continuum image obtained bv iMaurv et al.l (12011 


see their Figure 3). It is noteworthy that the local peaks of the HCO'*~ map a re coincident 
with or located close to massive dense clumps identified bv IPirogov et ahl (120131) through the 
1.2 mm continuum observations, although they mapped only a half of the area of our map. 
We found that the intense HCO’*' emission is seen mainly in two parts in the map: One is 
the arc-shaped structure in the south-west part, and the other is the filamentary structure 
in the north-east part. The central part of the cluster including IRS lA South ionizing the 
H II region is located in between the two parts. As seen in Figured] (b), YSOs are distributed 
around the dense gas traced in HCO’*', suggesting that the sources were formed therein. The 
HCO'''map also shows that there is less dense gas around IRS 5, though several Class O/I 
sources or starless cores have been detected there. 

We did not detect the CCS (JAr=43 — 32 ) emission in W40 at the present sensitivity. The 
CCS molecule is known to be niore abundant in an earlier sta ge of chemical evolution (e.g.. 


Suzuki et al.lll992l: IHirota et al.ll2009l: IShimoikura et al.ll2012h . It is interesting to note that 
the Serpens South cloud locat ed only ~ 20' to t he we st of W40 exhibits strong CCS emission 


(^mb ^ 2 K) as reported bv iNakamura et al.l (I2014J ). They estimated the column density 
and fractional abundance of CCS to be Y(CCS)~ 1 x 10^^ cm“^ and /(CCS)= 10“^° — 10“® 
in the Serpens South cloud. From the noise level of our CCS data, we estimate the upper 
limits of Y(CCS) and /(CCS) to be ~ 1.1 x 10^^ cm“^ and ~ 10“®'^^ at the intensity peak 
position of the HCO’*’ emission in W40, respectively (see Appendix]^. These values infer 
that Y(CCS) in W40 is likely to be an order of magnitude smaller than in the Serpens South 
cloud, which can be the cause of the non-detection of the emission line. The upper limit of 
/(CCS) is comparable to the fractional abundance in the Serpens South cloud, suggesting 
that W40 may be more evolved (or at a similar stage at most) in terms of the chemical 


"^Herschel SPIRE image which was obtained as part of Proposal ID:SDP_pandre_3. 
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reaction time. 

The velocity field of W40 shows a great complexity. Figure [2] displays the velocity 
channel maps of the and HCO’*'emission lines with a velocity interval of 0.4 km s“^. 

The emission is significantly detected in two velocity ranges: One is the range 3 < 

1/lsr ^ 6 km s“^ and the other is the range 9 < Id^sR ^ 12 km s“^. The emission is weak 
in the range 6 ;< Vlsr ^ 9 km s“^ where the HCO’*' emission is prominent. This indicates 
that there is a velocity component at Vlsr — 7 km s“^, and the optically thick emission 
line is heavily self-absorbed while the optically thinner HCO’*' emission appears as emission, 
which is obvious in the position-velocity (PV) diagrams shown in Figures |3] (a) and (b) 
together with the intensity-weighted mean velocity map of the HCO’*'emission in Figure [3] 
(c). Note that, in addition to the velocity component at ~ 7 km s“^, there are also distinct 
components at ~ 5 and ~ 10 km s“^ detected in HCO^ just around the Hll region. In the 
next subsection, we discuss the details of the velocity structures of the observed region. 


3.2. Four Velocity Components 

A careful inspection of the and HCO’*' data infers that the observed region consists 
of at least four velocity components at IAsr — 3, 5, 7, and 10 km s“^. As summarized in 
Table [H the ~ 7 km s“^ component represents the systemic velocity of the entire W40 region, 
and the ~ 5 and ~ 10 km s“^ components are likely to originate from dense gas located along 
the expanding shell of the H ll region. The dense gas corresponds to Clump 1 (~ 5 km s“^) 
and Clumps 5 and 6 (~ 10 km s“^) detected in HCO’*' that we will describe in Sections 13.51 
and 13.61 The ~ 3 km s“^ component is weak and shows up in a very limited region, and it 
might be unrelated to W40. As explained in the following, the four velocity components can 
be recognized in the ^^CO and HCO^ spectra. 

In FigurelU we show examples of the ^^CO and HCO^ spectra observed at some positions 
mainly around the YSOs together with the velocity channel maps with a coarser velocity 
interval (1.3 km s“^) than Figure 2. The four components can be seen in some of the 
spectra. For example, the spectra in panels (a)~(c) sampled in the west of the Hll region 
consist only of the ~ 5 km subcomponent, and those in panels (f) and (g) sampled in the 
south of the H ll region exhibit the ~ 10 km s“b component as well. Note that the ~ 5 and 
~ 10 km s“b components are detected both in b^CO and HCO^only around the Hll region. 
In addition, the HCO’*'spectra in panels (e) and (i) consist of the ~ 7 km subcomponent. 
The ~ 3 km s“b component is detected both in b^CQ and HCO’*' and is seen only around the 
positions shown in panels (j) and (k). 






Compared with the 2MASS image in Figure IU the ~ 5 and ~ 10 km components 
exhibit intense emission, and they are likely to delineate the expanding shell around the 
H II region. The ~ 7 km s“^ component is well detected in HCO’*' and is likely to represent 
the systemic velocity of the entire W40 region. As mentioned earlier, the absence of the 
intense emission for the ~ 7 km s“^ component indicates that there is a large amount 

of cold gas in the foreground of the Hll region, which should cause the heavy absorption 
in around this velocity. Actually, the emission line over the velocity range 

6 — 9 km s“^ have a brightness temperature of only Tmb — 2 K while the optically thinner 
HCOemission line is brighter (up to ~ 8 K) as seen in Figure H] (e). 


Grutcher 

1977; 

Piroerov et al. 

2013 

). Molecular observations done by 

PiroETov et al. 

(20L3) 


covers the arc-shaped structure just aro und the H ll regi on, and they detected the two velocity 
components at ~ 5 and ~ 7.5 km s“^. ICrutcherl (119771) made ^^CO (J = 1 — 0) observations 
with a lower angular resolution (HPBW= 2(6), and they detected the ^^CO emission at ~ 5 
and ~ 10 km s“^. The author also found an absorption line at ~ 7 km s“^ in the OH 
spectra obtained by a dditional o bserv ations with the NRAO 43 m telescope (HPBW= 18'). 
Based on these data, I Crutcher! (119771 ) concluded that W40 has a single component with a 
peak velocity at IAsr — 7 km s“^, which is extended in the wide velocity range 3 — 12 km 




s^" and is self-absorbed at around 7 km s“^ by colder gas in the foreground. We agree with 
their interpretation that W40 is surrounded by the ~ 7 km s“^ component. In addition, we 
found the existence of two other velocity components at ~ 5 and ~ 10 km s“^ around the 
H II region, probably originating from the dense gas located along the expanding shell, which 
are detected both in ^^CO (J = 3 — 2) and HCO’*' (J = 4 — 3) as seen in the PV diagram in 
Figure |3] as well as in the spectra in Figures 01 (b) and (f). 


3.3. Search for Outflows 

Looking at the ^^CO spectra in Figure 01 there are some cases in which high velocity 
wings are obviously seen in spite of the huge complexity in the velocity held. For example, 
the ^^CO spectrum in Figure 0] (a) exhibits a blue-shifted wing at Hlsr < 3 km s“^, and 
the one in Figure 0] (g) exhibits a red-shifted wing at I4sr > 12 km s“^. These wings 
are very likely to be tracing real high velocity gas, not mere Gaussian tails of the intense 
spectra of the ~ 5 and ~ 10 km s“^ components. As detailed in the Appendix |Bl most 
of the spectra in the observed region can be well htted by simple Gaussian functions 

except for the apparent high velocity wings as seen in Figure 0] (a) and (g). 

Because the high velocity wings are often seen in the ^^GO spectra along the dense 
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condensations traced in HCO’*' both aronnd the H ll region and the hlamentary strnctnre, 
we snggest that they are either the high velocity gas blown away by the expansion of the 
H II region or the ontflowing gas driven by YSOs forming in the dense condensations. At 
the moment, there is no way to distinguish these two effects reliably only with the present 
dataset, and both of them should be actually at play. Modeling the expected signature 
of instabilities in high velocity gas from H ll region expansion is beyond the scope of the 
current study. For the purpose of this study, we interpret apparent high velocity wings on 
the principal velocity components of as evidence of potential “candidate” outflows 

from YSOs and derive their parameters accordingly. Though a certain fraction of the high 
velocity wings can be due to the expansion of the H ll region, there are some high velocity 
wings found in the vicinity of YSOs (e.g., “B2” and “R2” in the following) including the 
one located away from the Hll region (e.g., “R4”) which are likely to be promising outflows 
from YSOs. 


Since the dense gas consists of the ~ 5 km s“^, ~ 7 km s“^, and ~ 10 km s“^ components, 
we searched for outflow candidates blowing out from these components. For the ~3 km 
s“^ component, we did not find any high velocity wings in the ^^CO spectra, and thus we 
disregarded this component. In order to carry out a systematic search for outflows, we set 
constant velocity ranges to hnd high velocity wings in the spectra. As illustrated in Fig- 
ureO we decided to use the velocity ranges —3.2 Y IAsr = 2.8 km s“^ and 11.7 Y Rlsr = 14.5 
km s“^ for the blue- and red-shifted high velocity wings of outflows, respectively. Basically, 
we can carry out a reliable search for outflow candidates in only in these velocity 

ranges, because the heavy absorption by the ~ 7 km s“^ component precludes the detection 
of the faint wing emission over the range 6 < IAsr ^ 9 km s“^. However, for a refer¬ 
ence, we will also investigate wing emission possibly associated with the ~ 5 and ~ 10 km 
components in the range 5.7 < Vt.sr ^ 8. 4 km s ~^ and 8.5 Y Vlsr = 9.2 km s“^, because. 


..-1 


in the vicinity of the H ll region, IZhu et al.l (120061 ) found a bipolar outflow associated with 
the ~ 5 km s“^ component, though we now realize that its red-shifted wing emission should 
be contaminated by the heavy absorption. 

Figure [6] shows the intensity distributions integrated over the velocity ranges 

mentioned in the above. The overall distributions of the high velocity gas in the figures 
show a great complexity similar to what has been reported for the Circinus cl oud where the 


complex outflow lobes are found to be a mixture of som e smaller outflows fe.g.. lDobashi et ah 


19981 : iBallv et al.lll999l: IShimoikura fc Dobashill201l[) . which is likely the case also for W40. 
In addition, the distributions in the figures are contaminated by some distinct velocity 
components (at Rlsr — 3 km s“^) unrelated to the outflows. 


The distributions on the sky are too complex to separate them reliably into single 
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outflow candidates only from the present dataset. In this paper, we thus decided to pick 
up the local peaks in the intensity map in Figure |6] if they are not due to the distinct 

velocity components but are very likely due to outflows showing a wing-like feature in the 
spectra. As a result, we found 13 local peaks as indicated in Figure El For simplicity, 
we will regard each of the 13 local peaks as the peak position of a single high velocity lobe 
and we will consider them to be candidates of outflows in W40. There are nine blue lobes 
and four red lobes. As summarized in Table [21 we call the blue lobes “B1”-“B9”, and the 
red lobes “R1”-“R4”. 

To measure some physical properties, we defined the surface area of the outflow lobes 
S'lobe at the 1.2 K km s“^ (corresponding to the ~ 9a noise level) contour level. In the case 
that two or more lobes are attached to each other at this threshold value, we divided them 
into single lobes at the valley in the contour maps of the intensity (Figure ED- We 

also attempted to find driving sources of the outflow lobes. It is however also difficult to 
identify the driving sources reliably, and thus we simply picked up YSOs located within S'lobe 
defined in the above as possible candidates for the driving sources. There are nine outflow 
candidates (Bl, B2, B4, B6, B7, B8, Rl, R2, and R4) with at least one candidate of the 
driving source, and there are four outflow candidates (B3, B5, B9, and R3) which could not 
be assigned to any known sources. 

We summarize the main features of the outflow lobes as follows. 


Blue lobes: Blue lobes are seen mainly around the H ll region. A candidat e for the driv¬ 


i ng so urce of Bl is 1.2 mm source, which is found and classified as “starless” bv iMaurv et ah 


(I 2 OIII) . A jet-like blue lobes B2 and B3 elongated in the East-West direction are clearly seen 
in Figure El B3 is the largest lobe located close to IRS 5. Multiple blue lobes B4-B8 are 
located to the south of the H ll region. B9 is a small lobe f ound in the filament structure. 
B1-B8 correspond to an outflow discovered by IZhu et ah! (120061. see their Figure 3) with 
lower angular resolution (80"), which is resolved into smaller multiple possible outflows in 
our map. 

Red lobes: Two red lobes Rl and R2 are detected in the vicinity of IRS lA South, 
suggesting that these outflow candidates possibly originate in t he cluster. C andidates of 


driving sources for these lobes are classified as Class I protostars flMallick et al.ll2013[) . The 
lobes R3 and R4 are found in the filament structure. 

Finally, in Figure [71(a) and (b), we compared the distributions of the outfl ow candidates 


with t he continuum-subtracted 2.12 /im H 2 narrowband image obtained bv [Mallick et ah 


fj2013[) to search for the possible counterparts of the CO outflow lobes. We also compared our 
map with the Spitzer IRAC 4.5 /im image in Figure[71(c) and (d) to see whether vibrationally 
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excited CO emission is detected around the outflow candidates. However, we could not 
identify clear counterparts corresponding to the outflow lobes in these images, though there 
are some near/mid-infrared nebulosity in the vicinity. 

According to the Herschelmap (the left panel of Figured] (b)), the Hll region is appar¬ 
ently expanding into the surrounding medium, and some of the identified lobes (i.e., B2, B3, 
Rl, and R2) are located in less dense regions at the edge of the Hll region. Furthermore, 
these lobes are associated with the ~ 5 and ~ 10 km s“^ components which are likely to rep¬ 
resent the expanding velocity of the shell (see Section jXSj) . For the lobes found just around 
the H II region, we cannot distinguish if the observed high velocity gas is due to real outflows 
from YSOs or is merely tracing the expanding shell. In addition to the above possibility, 
we should also note that some of the identified lobes might not originate from a single YSO 
(see Figure E]), but originate from the cluster in the Hll region. All of these ambiguities 
arise from the complex velocity filed in W40, especially the heavy absorption by the ~ 7 
km s“^ component. In this paper, we will regard the 13 blue- and red-shifted high velocity 
lobes to be candidates of outflows in this region and will estimate some outflow parameters 
in Section 13.41 


3.4. Properties of Outflow^ Candidates 


In order to estimate the masses of the outflow candidates, we assumed that the molecules 
are in the local thermal equilibrium (LTE) and also that the emission is optically thin 
at all wing velocities to calculate the column densities in the outflow lobes for a uniform 
excitation temperature of Tgx = 30 K. We then estimated the mass, momentum, and energy 
of the outflow lobes. Details of the calculations of the physical parameters are described in 
Appendix [C] and the results are summarized in Table |2j 


Zhu et al.l (120061) suggested a total mass of 0.48 Mq for the blue lobe of the outflow 


they found which corresponds to our B1 to B8, on the assumption of a distance of 600 pc 
and a mean molecular weight of 2.72 (cf., we assumed 500 pc and 2.4, respectively). The 
sum of the corresponding blue lobes in our Table [2] amounts only to ~ 0.11 Mq which would 
be rescaled to ~ 0.18 Mq for the same assumptions as theirs. The difference between their 
and our estimates is mainly due to the contamination of the distinct velocity component at 
WgR ~ 3 km s“^ (see Section which was not resolved in their observations, and also due 
to the fact that we defined the surface area of the outflow lobes at a rather high threshold 
value of 1.2 K km s“^ (see Section l33|) . 


As seen in the 8th column of Table [S] we found that the masses of the individual lobes 
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are in the range from ~ 0.001 Mq to ~ 0.03 Mq. The YSOs that can be the possible driv- 
i ng sources of thes e outflow candidates have bolometric luminosities Lboi — 4.2 — 6.6 Lq 


llManrv et al.l 1201 ll . see their Table 1). According to a statistical study of outflows by 
Wn et al.l (12004J ). outflows associated with low-mass YSOs with Lboi — 1 — 10 Lq have a 


mass range of 0.01 Mq to a few Mq. The masses of outflow candidates we found in W40 
are much smaller than the average mass of outflows associated with low-mass YSOs, and are 
close to the lower end in the plo ts of the outflow masses versus the bolometric luminosity of 
the driving sources reported by Wn et ah f 2004h (see their Figure 5). 


Here we estimate the possible errors in the derived parameters caused by our assump¬ 
tions. First, we assumed the constant excitation temperature Tgx = 30 K. An error arising 
from this assumption should be small, because the outflow masses are insensitive to Tgx, and 
they would fluctuate only by ~20% at most within the variation of Tgx in the mapped region 
(20 < Tgx < 65 K). Second, we assumed that wing components in the (J = 3 — 2) 

emission line are optical l y thin (r -C 1), but they might be moderately optically thick. For 
example, iDunham et al.l (1201411 investigated 17 outflows associated with low mass stars with 
the (J = 3 — 2) and (J = 3 — 2) emission lines, and found the optical depths of 

the (J = 3 — 2) emission line varies in the range r ~ 1 — 10 with an average value of 

~ 2.8. If we adopt their average value, the inverse of the escape probability in Equation flC3jl 
would be = 3.0. This infers that our optically thin assumption may cause an underes¬ 
timation in the mass estimate by a factor of a few. Finally, it is important to remark that 
the parameter Kvmg (i.e., the velocity range exhibiting the high velocity wings) in Equations 
(IC4p - (IC7p should be highly underestimated because of the heavy absorption in the velocity 
range 6 < IAsr ^ 9 km s“^, and hence our estimate for the mass, energy, and momentum 
of the outflows should be highly underestimated. 


To summarize our outflows survey, we found that 13 candidates of outflow lobes showing 
wing-like emission in the spectra. However, we should note that our estimates of the 

outflow parameters such as the mass, momentum, and energy should give only the lower 
limits to the actual values mainly because of the heavy absorption at ~ 7 km s“^ and the 



3.5. Clumps Detected in HCO 


The HCO’*' (J = 4 — 3) emi ssion line is a goo d tracer of cold dense gas for its high 


critical density > 10® cm ^(e.g. 


Chart et al.l l2000l) . As stated in Section [321 the radial 
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velocities of the HCO’*'emission can be classified to the three velocity components at ~ 5, 
7, and 10 km s“^ in Table [H and they can be mostly well separated at 14sr =6.0 and 9.0 
km s“^. We show the HCO’*' integrated intensity maps divided at these velocities in Figure 
IS In order to investigate physical properties of the dense gas detected in HCO+, we will 
dehne clumps based on these maps. We will regard a condensation enclosed by the lowest 
contours (1 K km s“^ ~ 6 a) having the peak integrated intensity greater than 2 K km s“^ in 
each map as one clump. There are hve condensations satisfying the criteria. However, a 
large condensation in the 6 ^ Wsr = 9 km s“^ range exhibits two slightly different velocity 
components by ~ 1 km s“^, and thus we divided them into two at the boundary of the 
two velocity components as labelled Clumps 3 and 4 in Figure | 8 l As a result, we identihed 
six clumps in the HCO’*' data. We show the HCO’*' spectra sampled at their peak positions 
in Figure m and also show the line parameters (i.e., the peak brightness temperature, line 
width, and peak velocity) obtained by htting the spectra with a single Gaussian function in 
Table H 

Based on the HCO^ data, we will derive molecular masses of the clumps on the assump¬ 
tion of the LTE using Equations flC2p and flC3p for X =HCO’'' (see Appendix [C]) . According 
to JPL CataloJ§ and SplatalogucI, the constants Bq, /r, and Eu are 44.5944 GHz, 3.888 
debye, and 29.7490 cm“^, respectively, for the HGO"*" (J = 4 — 3) emission line. We list these 
constants in Table |3l 


In order to derive the mass of the clumps, we estimated Tgx for HGO^ comparing our 
HGO"*" (J = 4 — 3) data with the HGO’*' (J = 1 — 0) data measured by Pirogov et ah (2013) 
toward some limited positions (see their Table 4). For these positions, we smoothed our 
HGO"*" (J = 4 — 3) spectra to the same angular resolution as the HGO"*" (J = 1 — 0) spectr a 
(~ 45"), and applied the non-LTE radiative transfer code RADEX (Ivan der Tak et al.ll2007n . 
Results of the analyse s infer Tpy >3 0 K to reproduce the hydrogen number density u(H 2 ) ~ 
1 X 10® cm“^ found bv iPirogov et al.l (120131) for these positions. This val ue of Tpy is consistent 
with the temperatures of dense gas around the H ll region measured by Pirogov et ah ( 20131) 
using NH 3 and dust continuum data (16—40 K, see their Table 7). In this paper, we therefore 
adopt Tex = 30 K throughout the observed region. 

If we assume that the emission line is optically thin (r 1), the HGO"*" column density 
can be derived by the following equation for Tex = 30 K, 


A^hco+ — 4.53 X 10 


11 


dv 


cm 


-2 


( 1 ) 


'Vra 


^http://spec.jpl.nasa.gov/ 
^http: //www.splatalogue.net 
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where the term f units of K km s“^. The velocity range for the integration 

frange IS the Same as shown in Figure El i.e., 3 ^ Vlsr = 6 km s“^, 6 ^ Vlsr = 9 km s“^, 
and 9 ^ Vlsr = 12 km s“^ for clumps at ~ 5 km s“^ (Clump 1), ~ 7 km s“^ (Clumps 2-4), 
and ~ 10 km s“^ (Clumps 5 and 6), respectively. 

Assuming the fract ional HCO’*'abundance ratio to be Ahco+/1Vh2 = 2.5 x 10“^° found 
by IPirogov et al.l (120131 ) through comparison of the HCO’*' (J = 1 — 0) emission line and 1.2 
mm continuum data toward our Clump 1, the total mass of the clumps can be estimated by 
the following equation, 


= 2.96 X 10-® D 


Z 


HCO+, 


mb 


v)dvdS Mq 


( 2 ) 


^ clump 


'v„ 


where the term J J Tj^J^'^^dndS' is in units of K km s“^ arcmin^ and D is the distance to the 
clumps in units of pc {D = 500 pc for W40). We dehned S'dump is the surface area of the 
clump at the 1 K km s“^ contour level. 

Using Equations ([1]) and ([2]), we estimated the total mass of the clumps. The results 
are presented in Table 01 As seen in the table, the identified clumps have a masses ranging 
from ~ 2 to 26 Mq with a mean value of ~ 12 M( 


©• 


Here, we should note that our mass estimates depend on the adopted Tex. The derived 
masses should increase for lower Tex, an d they would change by ~ 50 % at most if Tex varies 
in the range 16 — 40 K as measured by Pirogov et ah ( 20131) . 


Pirogov et ahl (1201311 suggested a total mass of 8.1 Mq (the sum of the masses of their 


clumps 2 and 3) for a part of our Clump 2 from the 1.2 mm continuum data. We estimate the 
mass contained in the same area measured by Pirogov et ah to be 9 Mq. Th e mass derived 


from our HCO"*" data is therefore in close agreement with that derived by iPirogov et ah 


(120131) . However, on a larger scale, iMaurv et al.l (120111) suggested a much larger mass than 
our estimate. They estimated the total mass of a larger region equivalent to our entire 
observed area to be 310 Mq through the Herschel column density map. They assumed a 
distance of 260 pc, and thus the mass should be rescaled to ~ 1150 Mq for our assumed 
distance (500 pc). We obtain only ~ 110 Mq as the total mass of the observed region, which 
is only 10% of their mass. This is probably due to the fact that the dust emission (detected 
by Herschel) traces all of the dust along the line-of-sight, while HCO’*'traces only very high 
de nsity regions becau se of its high critical density (> 10® cm“^). Moreover, as suggested 
bv iGirart et al.l (1200011 . HCO^ (J = 4 — 3) is not likely to be thermalized, and therefore the 
masses derived here are lower limits. 


Finally, we attempt to hnd the possible association between the identihed clumps and 
the outflow lobes found in the previous subsections, though it is difficult to establish their 
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definite association because of the enormous complexity in the velocity held for which we 
cannot identify even their driving sources reliably. Here, we regard that the clumps and the 
outhow lobes are associated if their surface areas overlap each other. All of the outhow lobes 
except for B1 and B3 have an overlap with one or more clumps. We will regard that B1 and 
B3 are associated with Clump 1, because weak HCO’*' emission is extending from this clump 
toward these outhow lobes. We also checked the systemic velocity 14ys for each of the lobes 
(see Appendix [O and the 6th column of Table [2]). For B7 and B8, these 14ys are close to the 
radial velocity of Clump 2 and we will regard these lobes are associated with Clump 2. The 
results summarized in the last column in Table HI 


Comparing the identihed clumps with the distributions of the YSOs in Figure [H It 
is interesting to note that some of the HCO+clumps (e.g., Clump 1 and Clump 4) are 
associated with Class O/I sources at their peak positions. Moreover, the clumps associated 
with the outhow lobes tend to have broader line widths (see Table ip, suggesting that the 
outhow lobes are inhuencing on the kinematics of the parent clumps. iNakamura fc Lil (1200711 
investigated theoretically the complex interplay between outhows induced by star formation 
and turbulence of the parent clumps, and showed that outhows should play a signihcant 
role in altering the star-forming environments. The energy and momentum of the outhows 
found in this study appear too small to give signihcant inhuence on the parent clumps, but 
again, we should note that our estimates for the outhow parameters may suher from a large 
uncertainty. 


3.6. Three-Dimensional Model of W40 


Here, we discuss the locations of the identihed clumps relative to the H ll region and 
propose an updated three-dimensional model of W40. Previous studies suggested that the 
H II r egion lies on the edge of a large molecular cloud in a support of the blister model 


Zeilik fc Lada 1978 


Valleel 1198711. More recently, Hersche l revealed the total dust 
Andre et ah 2010 : Konvves et ah 2010). We show the 250 /xm 


(e.g. 

distribution in W40 (e.g. 
image by Herschel and the three-colors {JHKs) composite image of 2MASS in Figure [9] 


overlaid with our HCO^ channel maps (~ 5 km s 


-1 


7 km s ^, and ~ 10 km s ^). 


Morphological coincidence between the molecular gas and the dust seen in the far-infrared 
emission [Herschel) and in the near-infrared absorption (2MASS) is obvious. It is noteworthy 
that Clumps 3 and 4 show up in the 2MASS image in the absorption, indicating that these 
clumps are located in the foreground of the H ll region. 


Among the four velocity components summarized in Table [H the ~ 7 km s ^ component 
is the primary component of molecular gas around W40. Our Clumps 2-4 exhibit this velocity 
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in HCO’*'. As seen in the PV diagram in Figure [3], there is a hea vy absorption fea t ure in 
the emission over the velocity range 6 < IAsr ^ 9 km s“^. iNakamura et ahl (120111 ) 

have shown that the same absorption feature in is clearly seen around t he Serpens 

South molecular cloud (~ 20' to the west of W40) at the same velocity range. iDzib et al. 


(I 2 OIOI) measured the distance to Serpens Main to be approximately 415 pc t hrough VLBI 


obser vations, which is much further than previously thought (e.g., 260 pc, iMaurv et ah 


20111). On the other hand, the distance of W40 has not been determined precisely varying 


in th e range 340-600 pc (e.g., iRodney fc ReipurthI l2008l: iKuhn et ahl I 2 OIOI: IShuping et al. 


20121) . and the W40 ar id Serpens South regio ns might be physically connected as suggested 


by some authors (e.g., iBontemps et ahl 120101) . A more likely explanation for the absorption 


is that a colder cloud, perhaps a part of the envelope surrounding both of the W40 and 
Serpens regions, lies in the foreground. We thus regard that Clumps 3 and 4 are a part of 
the larger cloud surrounding these regions, and that they are located in the foreground of 
W40. 


In the HCO^ data, the other velocity components at ~ 5 and ~ 10 km s ^ 
mostly around the Hll region. These components are very bright in (Tj 


- mb 


are seen 
> 40 K), 
to the 


indicating that they are probably interacting with the H ll region being locate d close 
ionization front. This picture is consistent with the blister model suggested bv I Valleel (119871) 


in which the Hll region is expanding away from the dense molecular co re located near the 
cloud surface. We will develop the picture suggested by IValleel (119871) . The ionized gas 
associated with the Hll region is expanding not only outward but also into the molecular 
cloud pushing away the surrounding dense gas. In this interpretation. Clumps 1, 2, and 
5 should be part of the expanding shell. It is noteworthy that these clumps have a radial 
velocity of ~ 5 km s“^, ~ 7.5 km s“^, and ~ 10 km s“^. This infers that they are the dense 
gas interacting with the Hll region being swept up by the expansion of the ionization front. 


On the basis of the above results and discussion, we suggest a three-dimensional model 
of the clump distributions in W40, which is summarized in the schematic illustration shown 
in Figure HO] In the following, we will discuss the locations of the individual clumps based 
on the HCO^ and ^^CO data and on their appearance in the 2MASS image (Figure E]). The 
molecular data provide us the information on the velocity and temperature of the clumps, 
and the scattered light seen in the 2MASS image often infers the geometry between the 
clumps and the H ll region. 

Clump 1: In the 2MASS image, this clump is located very close to the H ll region. The 
rim of the clump facing to the H ll region is bright, and the rest of the clump appears opaque 
(Figure El (b)). The ^^CO emission from this clump is very intense (~ 53 K, see Figure H] (c)) 
with a radial velocity of ~ 5 km s“^ which is blue-shifted by ~ 2 km s“^ compared to the 
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systemic velocity of the entire W40 region (~ 7 km s“^). We thns snggest that the clnmp is 
located on the near side of the expanding sphere of the H ll region being heated and swept 
np toward ns by the ionization front. 


Clump 2: This clnmp is the brightest in HCO’*' a mong all of the clu mps ( ~ 8 K, 
see Table 0]) with a high density of n(H 2 ) > 10® cm“^ flPirogov et al.N2013h . This clump 
overlaps with Clump 1 on the sky, but they can be clearly separated in velocity. We cannot 
identify Clump 2 on the 2MASS image, but its counterpart can be recognized in the Herschel 
250 yum image (see Figure ini(f)), suggesting that Clump 2 is located behind Clump 1. The 
radial velocity of Clump 2 is ~ 7.6 km s“^ which is consistent with the systemic velocity, and 
thus the c lump might be locat ed in the ambient background gas apart from the H ll region. 
However, iPirogov et al.l (120131 ) found a high kinetic temperature of 40 K for this clump 
based on the non-LTE modeling for the CS (J = 2 — 1) and (J = 5 — 4) emission lines (see 
their Table 7), suggesting that Clump 2 is probably in the vicinity of the Hll region. The 
above observational features would be naturally accounted for if we assume that the clump 
is interacting with the Hll region like in the case of Clump 1, being swept up toward the 
direction orthogonal to the line-of-sight, as illustrated in Figure [TUI fcb 


Clumps 3 and 4: These clumps are apart from the H ll region on the sky, and they 
can be easily identified in the 2MASS image as the opaque hlaments (Figure [9] (c)). The 
peak velocities of the HCO’*'emission line of Clump 3 is ~ 8 km s“^ and that of Clump 4 
is ~ 7 km s“^ consistent with the systemic velocity of the W40 region. These observational 
facts indicate that they are located in the foreground of the H ll region as explained earlier. 
However, we should note that Clumps 3 consists of two well defined peaks, and one of them in 
the south-west appears as a small clump with the faint scattered light in the 2MASS image, 
not as the opaque region. This part might be separated from the other peak with a different 
location along the line-of-sight, though they are connected by the weak HCO^ emission. 

Clump 5: This clump has a radial velocity of ~ 10 km s“^ which is red shifted by ~ 3 
km s“^from the systemic velocity. The clump is very bright in ^^CO (~ 50 K, see Figure 
S](f)), suggesting that the clump is interacting with the H ll region on the far side of the 
expanding sphere of the H ll region. The extent of Clump 5 coincides with the bright region 
in the 2MASS image around IRS lA South (Figure O (d)), and we cannot recognize the 
shape traced in HCO’*' in the image. 

Clump 6: This clump is located between the two local peaks of Clump 3 (at ~8 km s“^) 
on the sky, but it has a clearly different velocity (~ 9 km s“^). There is a faint nebulosity 
around Clump 6 in the 2MASS image. This clump is located rather far from the H ll region 
on the sky, and it exhibits ^^CO emission with a moderate brightness temperature of ~ 30 
K (Figured] (h)). We cannot decide the location of this clump with a conhdence, but we 
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suggest that it may be located rather in the back of the H ll region because of its appearance 
in the 2MASS image as well as of its slightly red-shifted radial velocity compared to the 
systemic velocity. 


As illustrated in Figure [101 Clumps 1, 2 , and 5 are likely to be located just around 
the Hll region. They are probably dense gas being compressed by the expansion of the 
HII region. It is interesting to note that 10 outflow lobes found in the data are likely 

to be associated with these clumps. This is suggestive of star formation triggered by the 
impact of expanding Hll regio n driven by nearby mass ive young stars, indicating that the 


sequential star formation (e.g.. lElmegreen &: Ladalll977D is taking place in this region. 


We also should note that three outflow lobes are likely to be associated with Clump 
4 which is away from the H ll region. It could be possible that the clump will be another 
site of cluster formation in W40. According to the study of 1 1 young clusters in the S247, 
S252, and BFS52 regions including the Monkey Head nebulae (IShimoikura et al.ll2013|) . four 
of the clusters actually consists of two or more smaller clusters with a separation of ~ 1 pc 
(e.g., see their Figure 11(c)). Such multiple cluster formation in a limited region may not 
be rare , because there are so me other young clusters showing multiplicity (e.g., the S235 
region, [Kirsanova et al.ll20141) . Though the relation between Clump 4 and the Hll region 


is unclear, the dust distributions on a larger scale indicate that there are two outstanding 
condensations of dust around Clump 4 and the Hll region (e.g., see Figure [T](b)), suggesting 
that the W40 region might be a site of sequential cluster formation. 


4. CONCLUSIONS 

We have carried out ^^CO (J = 3 — 2) and HCO+ (J = 4 — 3) observations around the 
H II region W40 with the ASTE 10 m telescope and identihed CO outflows and dense clumps 
in this region. We have also observed the region with molecular lines at 45 GHz such as SiO 
(J = 1 — 0, n = 0) and CCS (J7v=43 — 32 ) using the 45 m telescope at NRO. We summarize 
the main conclusions of this paper in the following points (l)-(4): 

(1) Intense ^^CO (up to Tmb — 50 K) and HCO+emission (up to ~ 8 K) lines were 
detected around the W40 H ll region. The velocity held in the region is highly complex, 
showing at least four distinct velocity components at 14 sr — 3, 5, 7, and 10 km s“^. The 
~ 7 km s“^ component represents the systemic velocity of the region which causes heavy 
absorption in the ^^CO spectra over the velocity range 6 < 14 sr ^ 9 km s“^. The ~ 5 and 
~ 10 km subcomponents exhibit high temperature (> 40 K), and observed mainly 

around the H ll region. The ~ 3 km s“b component is faint, and might be unrelated to W40. 
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(2) The SiO, CCS, HC3N, and HC5N emission lines were not detected at the present 
sensitivity. We estimated the npper limits of the colnmn density of CCS and its fractional 
abnndance of W40 to be ~ 1.1 x 10^^ cm“^ and ~ 10“®'^®, respectively. 

(3) Based on the ^^CO (J = 3 — 2) data, we made a survey for molecular outflows 
over the observed region. We searched for high velocity outflow lobes in the velocity range 
Vlsr = 2.8 km s“^ for the blue lobes and Wsr = 11-7 km s“^ for the red lobes. As a result, 
we identihed nine candidates for the blue lobes and four candidates for the red lobes, and 
estimated the lower limits of physical parameters such as the mass, momentum, and energy. 
Though some of the identihed lobes around the Hll regions (i.e., B2, B3, Rl, and R2) might 
not be real outhows but merely tracing the high velocity gas of the expanding shell of the 
H II region, it is difficult to distinguish the two phenomenon with the present dataset, because 
of the complex velocity held in W40. 

(4) Based on the HCO’*' (J = 4 — 3) data, we identihed six clumps and derived their 
physical properties such as the mass, radial velocities, and surface areas. We further discussed 
their locations relative to the H ll region, and suggested a three-dimensional model of their 
distributions. 
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gestions to improve this paper. We are also grateful to Dr. K. K. Mallick for his kindly 
providing us with the continuum-subtracted 2.12 /rm H 2 narrowband data. This work was 
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26610045, and 26287030) of Japan Society for the Promotion of Science (JSPS). The ASTE 
telescope is operated by National Astronomical Observatory of Japan (NAOJ) and the 45 
m radio telescope is operated by NRO, a branch of NAOJ. 
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A. Derivation of the Column Density and the Fractional Abundance of CCS 


We estimated the upper limit of A(CCS) using Equations flC3j) and flC4j) in Appendix 
El from the noise level of the obtained CCS data {a ~ 0.25 K km s ^ integrated over the 
velocity range 3 < Vlsr < 12 km s“^). We used constant s such as fi and Bq in the equations 
for the CCS molecule summarized by IShimoikura et al.l (120121 see their Table 4), and used 
Q of 24.302 for a typical excitation temperature of CCS (Tex = 5 K, e.g, iHirota et ah 


20091 1. To derive /(CCS), we smoothed the HCO’*'spectrum to the same angular resolution 
as CCS (57") at the HCO+intensity peak position (i.e., Clump 2 in Table 0]), and derived 


AfH?) from the HC O'*~ spectrum using the A(HCO’'' )-A(H 2 ) conversion factor reported by 


Pirogov et al.l (120131 ) in the same way described in Section IT^ We hnally estimated /(CCS) 


as A(CCS)/A(H 2 ). 


B. High Velocity Wings in the spectra 

In order to conhrm the existence of the high velocity wings seen in the ^^CO spectra, 
we htted the ~ 5 km s“^ and ~ 10 km s“^ components by simple Gaussian functions to see 
whether the high velocity wings are merely tracing outskirts of the intense ^^CO emission 
or not. We hrst masked the velocity range Wsr = 6 — 9 km s“^ where the ^^CO spectra 
are heavily contaminated by the absorption, and then htted each of the ~ 5 km s“^ and 
~ 10 km s“^ components with a Gaussian function with two components by applying a 
weight proportional to the second power of the observed brightness temperature at each 
velocity. We demonstrate some examples of the results in Figure [TTJ Panels (a), (b), and 
(c) of the hgure show the ^^GO spectra exhibiting high velocity wings around the ~ 5 km 
s“^ component, the ~ 10 km s“^ component, and no obvious high velocity wing, respectively, 
which are the same spectra shown in Figure ll](a), (g), and (h). As seen in Figure fTTl the 
^^GO spectra without apparent high velocity wings are nicely htted by the simple Gaussian 
functions, while those with the apparent wings cannot be htted well, indicating that the high 
velocity wings are not merely the outskirts of the intense ^^GO emission at ~ 5 km s“^ and 
~ 10 km s“^. 

We should also note that we performed the above Gaussian htting to all of the ^^GO spec¬ 
tra in the observed region, and subtracted the htted lines from the observed spectra to gen¬ 
erate a map of the high velocity wings. The resulting map (not shown) appears very similar 
to the blue- and red-shifted lobes shown in Figure [6l 
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C. Derivation of Outflows Properties 


In this appendix, we explain how we calculated the physical parameters of the outflow 
lobes. We flrst derived the excitation temperatures Tex of the observed region from the peak 
brightness temperature of the data at each observed position. In general, the observed 
brightness temperature can be expressed as 

= |J(T,J - J(Tb,)]|l - exp(-T)] (Cl) 

where J{T) = To/(e^°/^ — 1), Tbg = 2.7 K, and Tq = huo/k. Here, uq is the rest frequency of 
the (J = 3 — 2) emission line, k is the Boltzmann constant, and h is the Planck constant. 

The constant Tq is therefore 16.59 K for this emission line. We derived Tex from Equation 
fICip at each pixel in the map assuming that the line is optically thick (r ^ 1 ) at the 
velocities giving the maximum temperature to The resulting Tex varies depending on 

the positions with the maximum value 65.6 K and the average value ~ 30 K. We will assume 
Tex = 30 K to derive the physical parameters of the outflow lobes. 

In order to estimate the total masses of the outflow lobes, we assumed the local thermo¬ 
dynamic equilibrium (LTE). The column density of molecule X (e.g., X = or HCO’*'), 
Nx, can be derived from the observed spectra by the following equation, 

dNx = g T^^iv) dv (C2) 


where g is a. functi on of Tpx, the escape probability /5, and some molecular constants, and is 


expressed as fe.g.. iHirahara et al.lll992l) 


_ 3h Q 

^ Svrs e^o/Tex _ i _ j(Tbg) ' ^ ^ 

Here, Q is the partition function and is Q = kTe^/hB^ where Bq is the rotational constant 
of the molecule, /i is the dipole moment, Eu is the energy of the upper level, and Sij is the 
intrinsic line strength of the transition for i to j state. For the (J = 3 — 2) emission 

line, we used the values of 57.635 GHz, 0.11 debye, and 23.069 cm“^ for the constants Tq, /r, 
and Eu, respectively. We list these constants in Table [3l All of the constants were gathered 
from JPL Catalog and Splatalogue. (3 is related to the optical depth r as /3 = (1 — e~'^)/T. 
We assume that the high velocity wings of the (J = 3 — 2) emission line are optically 

thin (r 1), and therefore (3 = 1. 

For these values. Equation flC2p for X =CO leads to 


Nco = 4.36 X 10^^ 



T^^{v) dn 


cm 


(C4) 
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where the term f T^^dv is in units of K km s“^. The velocity range for the integration is 
—3.2 ^ Vwing ^ 2.8 km s“^ for the blue lobes and 11.7 ^ Kvmg ^ 14.5 km s“^ for the red 
lobes, as stated in Section 13.31 


The hydrogen molecular column density A^Ho is then calculat ed from Nco, assuming the 
fractional abundance iVco/A^H 2 = 1 x 10““^ (IFrerking et al.l 1198211 . The mass of the outflow 
lobes Mflow can be obtained by integrating the product amniVHj where mn is the hydrogen 
mass and a is the mean molecular weight assumed to be 2.4 corrected for the 20% helium 
abundance. 


The integration should be made over the surface area of the outflow lobes Fiobe as 


Mflow = 7.13 X 10-^ D 


T£°(v)dvdS Mr, 


0 


‘5'lobe 


(C5) 


where the term j f T^^dvdS is in units of K km s“^ arcmin^, and D is the distance to the 
outflow candidates in units of pc (D = 500 pc for W40). We defined S'lobe at the 1.2 K km 
s“^ (corresponding to the ~ 9a noise level) contour level as defined in Section [3731 


The total momentum Pflow and kinetic energy Pflow contained in the outflow lobes can 
be estimated by the following equations, 


Pflow = 7.13 X 10-^ P 




'S'lobe 


\4ysl duds' Mfz^kms 


-1 


0 ■ 


(C6) 


Pflow = 3.57 X 10"’'P 


- Vsysf dudS Mq km^ s 


“S'lobe ^-w 


2 ^-2 


(C7) 


where l/gys is the systemic velocity. In order to decide I^ys, we derived the HCO+ spectra 
averaged over the surface areas of each outflow lobe. We then htted the spectra with a single 
Gaussian function, and regarded the peak velocity of the Gaussian function as the systemic 
velocity I4ys. In Figure [121 we show the HGO'*'and ^^GO spectra averaged over the surface 
areas of the identihed high velocity lobes. The IIGO+ emission line of B3 is not detected, for 
which we assumed Ggys = 5 km s“^for the lobe. Using Equations flG5p - flG7p . we calculated 
^flow, Pflow, and Pflow for the individual outflow lobe. 


Next, we attempted to estimate the dynamical timescale, fd, of a outflow lobe. In 
general, td should be calculated using the maximum separation between the outflow lobes 
and the driving sources. For the outflow candidates found in W40, however, we cannot 
identify the definite driving sources and also cannot separate the individual outflow lobes 
reliably. We therefore estimated the dynamical timescale as p = Pflow/Tchar where Pflow is 
the radius of the lobes defined as Pflow = ^/SioheJ^ and Whar is the maximum difference 
between Uys and the velocities observed in the high velocity wings detected at > 3a (~ 0.5 
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K) noise level in the spectra observed at the intensity peak positions of the outflow 

lobes. We summarize the physical parameters of the 13 lobes derived in this appendix in 
Table O 
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Table 1. Velocity Components Observed in W40 


Vlsr 

(km s“l) 

Detection 
in HCO+ 

Detection 

in CO 

Comment 

~ 3 

Yes 

Yes 

Faint CO and HCO"*" emission which might be unrelated to W40. 

rw 5 

Yes 

Yes 

Intense CO emission (> 40 K) originating from 

the dense gas associated with the expanding shell of the H ii region. 

~ 7 

Yes 

Faint 

2 K) 

Systemic velocity component of the entire W40 region 
causing heavy absorption in the CO emission line. 

~ 10 

Yes 

Yes 

Intense CO emission (> 40 K) originating from 

the dense gas associated with the expanding shell of the H ii region. 






Table 2. Properties of the Candidates of Outflows 


Lobe 

Peak Position 
a(J2000) <5(J2000) 

m / ff'j 

‘5'lobe 

(arcmin^) 

/ / dvdS 

(K km s“^arcmin^) 

ksys 

(km s“^) 

^^har 
(km s“^) 

{Mq) 

^fiow 

(Mq km s“^) 

-^fiow 

(M0 km2 s-2) 

-^fiow 

(pc) 

id 

(10"* yr) 

B1 

18 31 11.0 

-2 02 45 

0.5 

1.0 

5.7 

5.5 

0.002 

0.01 

0.02 

0.06 

1.0 

B2 

18 31 11.7 

-2 03 45 

1.2 

3.8 

5.0 

6.7 

0.007 

0.03 

0.07 

0.09 

1.3 

B3 

18 31 19.0 

-2 03 05 

2.6 

12.6 

5.0 

9.3 

0.023 

0.08 

0.19 

0.13 

1.4 

B4 

18 31 11.0 

-2 06 45 

1.5 

4.2 

4.8 

5.5 

0.008 

0.03 

0.05 

0.10 

1.8 

B5 

18 31 11.7 

-2 08 15 

1.8 

7.5 

5.2 

5.3 

0.014 

0.05 

0.09 

0.11 

2.0 

B6 

18 31 15.7 

-2 07 25 

1.8 

9.9 

5.4 

5.0 

0.018 

0.06 

0.11 

0.11 

2.2 

B7 

18 31 19.7 

-2 06 45 

1.6 

5.8 

7.4 

6.4 

0.011 

0.06 

0.18 

0.11 

1.6 

B8 

18 31 22.3 

-2 07 25 

2.0 

14.3 

7.3 

8.3 

0.026 

0.14 

0.41 

0.12 

1.4 

B9 

18 31 53.7 

-2 02 25 

0.7 

1.2 

6.6 

6.1 

0.002 

0.01 

0.03 

0.07 

1.1 

R1 

18 31 26.3 

-2 05 25 

1.5 

4.9 

9.8 

4.2 

0.009 

0.02 

0.03 

0.10 

2.3 

R2 

18 31 31.7 

-2 05 25 

1.3 

3.8 

10.1 

4.3 

0.007 

0.02 

0.02 

0.10 

2.1 

R3 

18 31 40.3 

-2 02 05 

0.5 

0.8 

7.1 

6.3 

0.002 

0.01 

0.02 

0.06 

0.9 

R4 

18 31 48.3 

-2 02 05 

0.3 

0.5 

6.8 

6.6 

0.001 

0.01 

0.01 

0.05 

0.7 


Note. — The velocity range Vwing used in Equations llC4ll - llC7ll to calculate f is set to —3.2 ^ Vlsr = 2.8 km s ^ for the blue lobes (B1-B9), and 

11.7 ^ I^lsr ^ 14.5 km for the red lobes (R1-R4). 








Table 3. Molecular Constants 


Molecule 

Transition 

Sij 

Bo 

(GHz) 

(Debye) 

Eu 

(cm~^) 


II 

CO 

1 

to 

3.0 

57.6350 

0.110 

23.0690 

HCO+ 

J = 4-3 

4.0 

44.5944 

3.888 

29.7490 


Note. — All of the constants are taken from JPL Catalog and 
Splatalogue. 
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Table 4. Properties of the Clumps 


clump 

Peak Position 
a(J2000) 5(J2000) 

m / //^ 

Tlnb 

(K) 

UsR 

(km 

AU 

(km s“^) 

lV(HCO+) 

(cm-2) 

‘5'clump 

(arcmin^) 

//T^bd*;dS 
(K km s“^arcmin^) 

-^clump 

(Mq) 

Associated 

Outflow Lobes 

1 

18 31 19.0 -2 07 05 

3.5 

5.7 

3.2 

5.5x1012 

8.4 

19.3 

14.3 B1,B2,B4,B5,B6 
(B3) 

2 

18 31 21.0 -2 06 25 

8.1 

7.6 

2.3 

l.lxlO^® 

5.3 

24.0 

17.8 

B7,B8 

3 

18 31 33.7 -2 04 15 

4.5 

8.0 

1.7 

3.8x1012 

3.4 

10.2 

7.6 

(R3) 

4 

18 31 46.3 -2 04 25 

6.4 

6.9 

1.5 

5.1x1012 

14.4 

35.4 

26.2 

B9,R4 

(R3) 

5 

18 31 26.3 -2 05 45 

1.6 

10.3 

2.3 

1.7x1012 

2.8 

4.6 

3.4 

R1,R2 

6 

18 31 36.3 -2 03 55 

3.4 

9.1 

1.7 

2.9x1012 

1.5 

2.7 

2.0 



Note. — The mass of Clump No. 6 may be underestimated by ~ 50%, because of the fixed velocity range used to define the 
clumps (9 ^ Vlsr = 12 km s“^). The other half of the mass is included in the Clump 3. The outflow lobe B3 is outside of 
the surface area of Clump 1, but we regarded them to be associated with each other, because they are close to each other on 
the sky. Another outflow lobe R3 is located right in the middle of Clumps 3 and 4, and we cannot judge which clump R3 is 
associated with. 
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(a) i2CO(J=3-2) + DSS2-Red image 
O ,- 

map resolution 0-5 pc 


scale for'2CO(J=3-2) 

15.0 77.5 140.0 (Kkm/s) 



18 ^ 31 - 40 ^ 20 * 

Right Ascension (2000) 


(b) HCO+(J=4-3) + Herschel/SP\RE 250 |im image 




00 ' 

02 ' 

04' 

06' 

-2°08' 


18’32'40' 00' 3r20‘ 30'40‘ 

Right Ascension (2000) 


18’3r40^ 20“ 

Right Ascension (2000) 


00 ' - 


• 2 = 20 ’ 


O starless cores from Maury et al. (2011) 

+ Class 0 or Class 0/1 sources from Maury et al. (2011) 
O Class I sources from Maury et al. (2011) 

O Class I sources from Mallick et al. (2013) 


Fig. 1.— Integrated intensity maps of the (a) (J = 3 — 2) emission line and (b) 

HCO'*' (J = 4 — 3) emission line around W40 obtained with the ASTE telescope. Velocity 
ranges used for the integration are —3.0 A 14 sr = 15.0 km s“^ and 3.0 ^ IAsr = 11-4 
km s“^ for the and HCO^ maps, respectively. The lowest contours and the contour 

intervals for the map in panel (a) are 10 K km s“^, and those in panel (b) are 1 K km s“^. 
The red color scale in panel (a) is the DSS2-Red image, and the gray scale in panel (b) is 
the Hersch el SPIRE 250 um im age. Positions of the high mass IRS sources (IRS lA Sourth 
and IRS 5, Shuping et al. 2012h are indicated by star marks. Squares, plu s signs , and circles 
denote the starless cores a nd Cl ass O/I sources classihed bv iMaurv et ahl (120111 ) (colored in 
green) and iMallick et al.l (120131) (colored in blue). The map resolution (31") and a linear 
scale of 0.5 pc (at 11=500 pc) are shown at the top of panel (a). The area enclosed by the 
thin lines in panel (a) and right side of panel (b) denotes the observed region. 
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Fig. 2.— Channel maps of the (J = 3 —2) (gray) and HCO'*' (J = 4 —3) (red) emission 
lines in the velocity range 0.2 ^ Clsr = 13.0 km s“^ with the 0.4 km s“^step. The velocity 
range in units of km s“^ used for the integration is indicated at the top right corner of each 
panel. The lowest contours and the contour intervals are 3 K km s“^ for the maps, and 
are 0.3 K km s“^ for the HCO"*" maps. 
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(a) 


(b) 



O 

o 

o 

c 
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(c) HCO+{J=4-3) 


’^CO (J=3-2) 
HCO*(J=4-3) 


IS'-aiMO’ 20^ 

Right Ascension (2000) 


Fig. 3.— (a) Position-velocity diagrams of the (J = 3 — 2) (black) and HCO+ (J = 
4 — 3) (red) emission lines taken along the cut labeled (a) in panel (c). (b) Same as panel 
(a), but for the other cut in panel (c). Contours in the position-velocity diagrams are drawn 
at 0.4, 0.8, 1, 2, 4, 8, 10, 20, 30, and 40 K for the emission line, and 0.2, 0.4, 1.0, 
1.2, 1.4, 2.0, 2.2, and 2.4 K for the HCO^ emission line, (c) Mean velocity map of the 
HCOemission line over the velocity range 3.0 ^ Clsr = 11-5 km s“^. Contours represent 
the HCO’*'integrated intensity map same as in Figure [T](b). 
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Fig. 4.— Channel maps of the (black contonrs) and HCO’*' (red contonrs) emission 

lines in the velocity range 1.3 ^ 14 lsr = 13.0 with the 1.3 km s“^ step. The velocity range 
in nnits of km s“^ used for the integration is indicated in the top right corner of each panel. 
The lowest contours and the contour intervals are 5 K km s“^ for the maps, and 1 K 

km s“^for the HCO^ maps. The 2MASS three-band composite image {JHKs) is overlaid 
in the last panel. The symbols for the YSOs are the same as in Figure [H Representative 
(black lines) and HCO+ (red lines) spectra taken at some positions are also shown in 
panels (a) to (k). The HCO+spectra are scaled up by a factor 2. The vertical broken lines 
indicate the typical velocities of the three components at Vlsr —5 km s“^, 7 km s“^, and 
10 km s“^ in Table [B 
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Fig. 5.— Schematic illustration to search for the high velocity wings of outflows in the 
spectra. In the W40 region, there are three main velocity components at F^sr — 5, 
7, 10 km s“^ (Table [T]) , but the heavy absorption by the ~ 7 km s“^ component precludes 
the detection of the high velocity wings over the velocity range 6 < F^sr 9 km s“^. 
We therefore searched for the emission in the velocity ranges —3.2 ^ Vlsr = 2.8 km 

s“^ and 11.7 ^ Vlsr = 14.5 km s“^ to find candidates of the blue- and red-shifted wings of the 
outflows, respectively. The intensity integrated over these velocity ranges are shown 

by the blue and red solid contours in Figure |6l Outflow candidates originating from the ~ 5 
km s“^ and ~ 10 km s“^ components may show symmetric wing components, and so, for a 
reference, we further investigated possible wing components contained in the velocity ranges 
5.7 ^ I4 sr = 8-4 km s“^ and 8.5 ^ Vlsr = 9.2 km s“^. The intensities integrated over 
these velocity ranges are shown by the red and blue broken contours in Figure |6l 












36 


(a) 



Right Ascension (2000) 


(b) 



Fig. 6.— (a) Locations of blue lobes of the outflow candidates which are labeled “Bl”- 
“B9”. White plus signs denote the peak positions of the lobes. Blue solid and red broken 
contours represent the intensity integrated over the velocity range —3.2 ^ Id^sR = 2.8 

km s“^ and 5.7 ^ Wsr = 8-4 km s“^, respectively. The former contours are to search for 
the blue lobes of outflows, and the latter contours are made for a reference (see text and 
Figure |5]). The lowest contour and the contour interval are 1.2 K km s“^, and symbols for 
the YSOs are the same as in Figured! Clumps indicated by thin solid circles with arrows are 
due to the distinct velocity component (at ~ 3 km s“^) unrelated to the lobes of the outflow 
candidates, (b) Same as (a), but for red lobes of the outflow candidates. Identihed lobes 
are labeled “R1”-“R4”. Red solid and blue broken contours represent the intensity 

integrated over the velocity range 11.7 Y Wsr = 14.5 km s“^ and 8.5 Y Wsr = 9-2 km s“^, 
respectively. The former contours are to search for the red lobes of outflows, and the latter 
contours are made for a reference (see text and Figure [5|). 
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grayscale: Continuum-subtracted 2.12^m Hg narrowband image 

(a) (b) 



18’3r40* 20‘ 

Right Ascension (2000) 


18^3r40- 20‘ 

Right Ascension (2000) 


grayscale: Spitzer IRAC 4.5Mm image 

(C) (d) 



18’3r40' 20‘ 18'3r40‘ 20' 

Right Ascension (2000) Right Ascension (2000) 


Fig. 7.— Comparison of the CO outflow candidates (contour) and other emission lines (gray 
scale), (a) Locations of blue lobes of the outflow candidates (same as in Figure EKa)) overlaid 


with t he continuum-subtracted 2.12 /im H 2 narrowband image obtained by iMallick et al. 
(120131) . (b) Same as (a), but for red lobes (same as in Figure [H](b)). (c) Locations of blue 
lobes of the outflow candidates overlaid with the Spitzer IRAC 4.5 /rm image, (d) Same as 
(c), but for red lobes. The symbols for the YSOs are the same as in Figure [1] 
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Fig. 8.— Top panel shows the total integrated intensity map of HCO’*' same as Fignre^b) 
overlaid with the candidate outflow lobes. The other three panels below the top panel display 
the channel maps for the ~ 5 km s“^, ~ 7 km s“^, and ~ 10 km s“^ components integrated 
over the range 3 ^ Id^sR = 6 km s“^, 6 ^ hisR = 9 km s“^, and 9 ^ Id^sR = 12 km s“^, 
respectively. Six identihed clumps are indicated in these panels. The lowest contours and 
the contour intervals in the three panels are 1 K km s“^. Peak positions of the six clumps 
are indicated by red plus signs. The symbols for the YSOs are the same as in Figured] The 
HCO'*' spectra sampled at the peak positions are also shown. 
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(a) 2MASS 3 band composite image 
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(c) 2MASS + ~ 7 km e’ component 
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(e) Herschel/SP\HE 250 pm image 
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Fig. 9.— Distributions of dense gas traced in HCO"'" and dust revealed by 2MASS and 
Herschel. Panel (a) displays the three-colors {JHKs) composite image of 2MASS. Panels 
(b)-(d) display the same 2MASS image, but overlaid with the HCO'*' intensity maps (con¬ 
tours) integrated over the velocity range 3 ^ IAsr = 6 km s"'^ (blue contours), 6 ^ IAsr = 9 
km s”'^ (black contours), and 9 ^ IAsr = 12 km s“^ (red contours), respectively. Panel (e) 
displays the Herschel 250 pm image. Panel (d) displays the same Herschel image, but over¬ 
laid with the HCO’*' intensity maps. Blue plus signs denote the locations of IRS lA South 
and IRS 5. 











40 



IRS 1A South 


(c) 



Right Ascension 


line of sight 



<- 

Right Ascension 


(b) 


I 



'Symbols* 





ambient gas / N 
(-7 km s’) 1 J 

Hll region 

cluster 


outflow 

QoP 

clump 



■ 


* 

5 km s ‘ 

- 7 km s ’ 

10 km s ’ 


Fig. 10.— Schematic illustration of the proposed model for the three-dimensional geometry 
of the W40 region. As indicated in the box drawn by the broken line, panel (a) shows the 
view as is observed on the sky (e.g., see Figure [9]^b)-(d)), and panel (b) shows the view 
observed from the declination axis orthogonal to the line-of-sight. Panel (c) shows the close- 
up view around the H ll region observed from the right ascension axis (Clumps 3 and 4 are 
not drawn). In the model, Clump 1 (IAsr — 5 km s“^) is located on the near side of the 
expanding sphere of the H ll region, and Clumps 5 and 6 (~ 10 km s“^) are located on the 
far side of the same sphere. Clumps 3 and 4 (~ 7 km s“^) are located in the foreground of 
the H II region. 
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Fig. 11.— Examples of the Gaussian functions (blue lines) best fitting the observed 
emission lines (black lines) exhibiting (a) a blue-shifted high velocity wing associ¬ 
ated with the ~ 5 km s“^ component, (b) a red-shifted high velocity wing associated with 
the ~ 10 km s“^ component, and (c) no apparent high velocity wing. Each of the ~ 5 and 
~ 10 km subcomponents is htted with two Gaussian functions. The spectra in panels (a), 
(b), and (c) are the ^^00 spectra shown in panels (a), (g), and (h) of Figure 01 respectively. 












-42- 











Spectrum 


HCO’(J=4-3)-‘5 


Fig. 12.— The (J = 3 — 2) (black) and HCO+ (J = 4 — 3) (red) spectra averaged 

over the surface areas of the identihed outflow lobes. The HCO'*’ spectra are scaled up by a 
factor of 5. 








































































